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We report on a study of B mesons on Nf = 2 full QCD configurations using an RG-improved gauge action, 
NRQCD heavy quark action and tadpole-improved clover light quark action. Results on the heavy-light spectrum 
and the decay constants from 16 3 x 32 lattices at a -1 ~ 1.5 GeV are presented, and compared with quenched 
results obtained with the same action combination at matching lattice spacings. 
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1. Introduction 

The decay constant fg is being studied exten- 
sively on the lattice because of its importance for 
the determination of CKM matrix elements. The 
spectrum of excited B mesons and b baryons is 
being measured in present experiments, whereas 
there exist only few lattice results on this subject. 

In this article we report on our study of B 
mesons in two-flavor full QCD employing the 
NRQCD action for heavy quark and a tadpole- 
improved clover action for light quark. The dy- 
namical configurations have been generated using 
the same light quark action and an RG-improved 
gauge action with a plaquette and a rectangu- 
lar term. Details on our full QCD configurations 
can be found in Refs. JT|. A parallel study of B 
mesons using the clover action for heavy quark is 
presented in Ref. 

2. Simulation Details 

We present results for two sets of dynamical lat- 
tices corresponding to the heaviest and the light- 



Table 1 

Parameters of lattices. The statistics for the dy- 
namical lattices has been increased since Lat- 
tice'99. The scale is fixed by *Jd = A27 MeV 
(for each sea quark for dynamical configurations) 



^sca 


0.1375 


0.1410 


oo 


mps/my 
a- 1 [GeV] 
#conf. 


0.8048(9) 
0.937(6) 

648 


0.586(3) 

1.127(10) 

490 


0.919(7) 
195 
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est sea quark in our configuration set at f3 = 1.95. 
The results are compared to those from quenched 
lattices generated with the same RG-improved 
gauge action at f3 = 2.187, the lattice spacing 
from the string tension matched to the dynami- 
cal lattice with n sea = 0.1375. Some details on 
these runs are given in Table @. 

We take 5 k values for the light valence quark 
corresponding to mps/my ~ 0.8 — 0.5. The 
strange quark mass m s is fixed using the K and 
the <f> meson. Our results for the B s meson are 
obtained with m s from the K, and the </> is used 
to estimate the systematic error. 
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Table 2 

Results for decay constants. Errors given in this 
table are statistical (including the statistical un- 
certainty in Mb), and, where applicable, the un- 
certainty in fixing the strange quark mass. Other 
systematic errors are discussed in the text. 



^sca 


/b [MeV] 


/sjMeV] 


Jb s 1 }b 


OO 


193(4) 


221(4)(+7) 


1.147(10)(35) 


0.f375 


216(4) 


250(4)(+8) 


1.157(9)(+35) 


0.1410 


215(6) 


251(6)(+6) 


1.166(14) (+31) 



For the heavy quarks, we use NRQCD at 
0(1 /M) with a symmetric evolution equation as 
defined in ||. We employ 5 bare heavy quark 
masses, covering a range of roughly 2.5 — 4.5 GeV. 

The heavy-light meson mass M is determined 
from the difference of the meson energy at finite 
momentum and at rest, assuming the dispersion 
relation, E(p) - E(0) = ^p 2 + M 2 - M. As a 
consistency check, we use both the Bd and the 
B s meson to determine the b quark mass. 

In our calculation of decay constants, the 
heavy- light current is corrected through O(afM). 
The mixing coefficients between the lattice oper- 
ators Q contributing at this order to the time 
component of the axial vector current J4, and 
the matching factor to the continuum current has 
been calculated J3| in one-loop perturbation the- 
ory, 



Ji 



(l+ap Q )4J at + (l + a Pl )4 
7(2) 



(i) 

lat 



(i) 



For the RG-improved gluon action, ay has 
not been calculated, and we use a tadpole- 
improved one-loop expression for the MS cou- 
pling, a jjgi l l a )- 

3. Decay Constants 

Our preliminary results for fg s and /s s / fs 
are given in Table |[ along with the statistical 
error and, where applicable, the uncertainty in 
the determination of m s . Additional systematic 
errors are estimated as follows: 0(a 2 ) correc- 
tions, taken to be a 2 x O(l), are 5%. A previous 
NRQCD calculation using the plaquette gluon ac- 
tion at a" 1 - 1 GeV finds the tree level 0(1/M 2 ) 
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Figure 1. $ ee (a s (M)/a s (M B ))( 2 / /3 «)/ % /M 
(top), and one-loop corrections to fVM (bot- 
tom) as a function of the inverse pseudoscalar 
meson mass. In the upper plot, squares stand for 
K SC a = 0.1375, diamonds for /c sca = 0.1410 and 
fancy squares for quenched. In the lower plot, 

circles denote ap J$ at / ' J 4 , squares, otp\J^ J at / J4, 
(2) 

and diamonds, ap\ J\ lat / 'J4. 



corrections to be ~ 2% we estimate our er- 
ror from the truncation of the 1 /M expansion to 
be ~ 4%. The leading discretization effects from 
the light quarks and gluons of 0(aah.QCD) and 
0(a 2 Ag CD ) are 5%. Added in quadrature, these 
estimates give 7%. 

Our two-flavor results for fg and fs 3 given in 
Table || show a 10% increase compared to the 
quenched values (see also Fig. ||). We do not re- 
solve any sea quark mass dependence. The de- 
pendence on the value of a s is weaker than for 
the plaquette gauge action, and the difference be- 
tween renormalized and bare decay constants is 
only about 5%. 

In Fig. [| we show the one-loop corrections to 
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the current J4 as a function of the heavy-light me- 
son mass. In the B region, 1 /M ~ 0.2, we find the 
correction to J^} at to be very small and the two 

other terms to contribute about the same amount. 
(2) 

The J4 f t contribution also contains a discretiza- 
tion correction to the current first pointed out 
in H . We note that this discretization correction 
is considerably smaller for the RG gauge action 
than for the plaquette gauge action ||. 

For fsj fs, we cannot resolve a difference be- 
tween the three lattices. 

In a parallel study of B mesons using clover 
heavy quarks Q, we have obtained fs and fs s 
taking the chiral limit for sea quark at /3 = 
1.8, 1.95 and 2.1. The results from that study at 
(3 — 1.95 agree within the estimated errors with 
the present results from NRQCD. 

4. Spectrum 

In Fig. ||, we give our results for several B 
splittings from the lattices with nt = and 
rif = 2, n sea — 0.1375. The top part of the figure 
shows the B* — B splitting. At present, we can- 
not resolve any unquenching effects. For quarko- 
nia, on the same lattices, the hyperfine splitting 
is found to increase from the quenched value only 
by a few MeV §. We find the B* - B splitting 
to be ~ 30% smaller than the experimental value. 
Possible sources of systematic error are the finite- 
ness of the sea quark mass, the 0(a) correction 
to the coefficient of the a ■ B operator, and higher 
order relativistic corrections. 

In the middle part of Figure ||, we show re- 
sults for the B% — B splitting, and in the lower 
part, the spin-averaged At, — B splitting. We do 
not find significant unquenching effects. However, 
for definite conclusions, we need to study several 
sea quark masses and lattice spacings, which is in 
progress. 
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TM are supported by the Research for the Fu- 
ture Program of JSPS. 
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Figure 2. Meson and baryon splittings in 
MeV. Circles denote results from CP-PACS at 
a _1 ( v / o : ) ~ 0.9 GeV. Diamonds stand for results 
from (quenched) and || (n/ = 2). Only statis- 
tical errors are shown. The solid line denotes the 
experimental value, the dashed lines, its error. 
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